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Abstract. We study the perspectives of resonant and nonresonant charmed meson production in p + A reac- 
tions within the Multiple Scattering Monte Carlo (MSMC) approach. We calculate the production of the resonances 
!f - (3770), 1^(4040) and ^(4160) on various nuclei, their propagation and decay to D, D, D* , D* , D s , D s in the 
medium and vacuum, respectively. The modifications of the open charm vector mesons in the nuclear medium are 
found to be rather moderate or even small such that dilepton spectroscopy will require an invariant mass resolution of 
a few MeV. Furthermore, the elastic and inelastic interactions of the open charm mesons in the medium are taken into 
account, which can be related to (u, d)-, s- or c-quark exchange with nucleons. It is found that by studying the D/D 
ratio for low momenta in the laboratory (< 2 — 2.2 GeV/c) as a function of target mass A stringent constraints on the 
c-quark exchange cross section can be obtained. On the other hand, the ratios Dj/TJj" as well as D/DJ and D/Df 
at low momenta as a function of A will permit to fix independently the strength of the s-quark exchange reaction in 
Z?7 N scattering. 
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1 Introduction 



In the last decade the dynamics of charm quark degrees of 
freedom have gained sizeable interest especially in the con- 
text of a phase transition to the quark-gluon plasma (QGP) 
where charmed meson states should no longer be formed due 
to color screening |Jj],§,^,^^]. However, the suppression of the 
cc vector mesons in the high density phase of nucleus-nucleus 
collisions might also be attributed to inelastic comover scat- 
tering (cf. and Refs. therein) if the corresponding 
charmonium-hadron cross sections are in the order of a few 
mb \ 10, HI]. Also the J/<PN cross section is not known suffi- 
ciently well since the photoproduction data suggest a value 
of 3-4 mb WW, while the charmonium absorption on nucleons 
(at high relative momentum) in p + A and A + A reactions 
is conventionally fitted by 6-7 mb [[^[l4|] which might relate 
to the cross section of a pre-resonance charmonium state with 
nucleons before the actual hadronic state is formed. Present es- 
timates for the J/\P formation times range from 0.3 - 0.5 fm/c 
1 13], but additional data from p induced reactions on nuclei, 
where charmonium can be formed on resonance with moder- 
ate momenta relative to the tar get nucleons, should help to pin 
down the present uncertainties []15j], Such studies might be per- 



formed experimentally at the HESR, which is presently pro- 
posed as a future facility at GSI [|l6j]. 

In a previous study we have explored the perspectives of 
charmed meson - nucleon scattering in the pd reaction [ |l7| l 
where the charmonium is produced on resonance with the pro- 
ton in the deuteron and may scatter with the spectator neutron. 
Another possibility is also the reaction ir + d — ► J/ippp as dis- 
cussed by Brodsky and Miller in Ref. fll8|]. Here we investi- 
gate pA reactions and study the production of the resonances 
!?(3770), ^(4040) and <?(4160)_in various nuclei, their propa- 
gation and decay to D, D, D* , D* 1 D s , D s in the medium and 
vacuum, respectively. Furthermore, the elastic and inelastic in- 
teractions of the open charm mesons in the medium are fol- 
lowed up in the Multiple Scattering Monte Carlo (MSMC) ap- 
proach to study the dominant medium effects as a function of 
the target mass A. 

Our work is organized as follows: In Section 2 we will 
briefly describe the MSMC approach and evaluate the forma- 
tion cross sections for the resonances !? r (3770), tf'(4040) and 
^(4160) on protons and nuclei. The fractional decay of these 
resonances to open charm mesons in the medium and vacuum, 
respectively, is calculated in Section 3 whereas the dynamics of 
the open charm mesons is described and studied in Section 4. 
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Section 5 concludes this paper with a summary and discussion 
of open problems. 



2 Resonance production and decay in pp 
and pA reactions 

We here examine the possibility to measure the in-medium life 
time (or total width) of the resonances V = (#(3770), #(4040) 
#(4160)) produced in pA reactions. Before doing so one has 
to determine first the resonance properties in vacuum, i.e. their 
production in pp reactions. To this aim we describe the vector 
meson spectral function by a Breit-Wigner distribution 



A V (M 2 ) = N v 



Mr? ot {M) 



(M 2 



M 2 



$tn v (M)) 2 + M 2 r? ot (M) ' 
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where My denotes the vacuum pole mass, SRTTy is the vector 
meson self energy in the medium - which vanished in the vac- 
uum - and Ny ~ l/V (for small width r to t) is a normalization 
factor that ensures / dM 2 A{M 2 ) = 1. The total width f tot is 
separated into vacuum and medium contributions as 



where 



r tot (M) = r vac (M) + r coll (M), 



r vac {M) = Y j r c {M) 



(2) 



(3) 



with r c (M) denoting the partial width to the decay channels 
c = DD, DD* ,D*D, D*D*,D S D S ,D* S D S ,D S D*,D* S D*, re- 
spectively. The in-medium collisional width is determined from 
the imaginary part of the forward vector-meson nucleon scat- 
tering amplitude as 



r, 



coll 



47T 
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where p c is the meson momentum in the rest frame of the res- 
onance. In case of the #(3770) only the lowest channel to 
DD contributes such that the matrix element in (Q) can be 
fixed by the total width at half maximum which we take as 
r F WHM = 25 MeV [glp. For the #(4040) and #(4160) we 
adopt rFWHM = 50 MeV and FpwHM = 80 MeV, respec- 
tively [|TJ. The branching ratios for the #(4040) show a strong 
anomaly favoring the decays to channels with vector mesons 
D*. The PDG [fl]] quotes r(D Q D°)/r(D*°D° + c.c.) = 
0.05 ± 0.03 and r(D* D*°)/r(D* a D a + c.c.) = 32 ± 12, 
where the phase space factor ~ p\ already has been divided 
out. The origin for such a large variation in the matrix elements 
is not known so far. We use the actual numbers (within the 
range of uncertainty) 



where pa denotes the nuclear density. Furthermore, Sy/y (0) is 
determined via the optical theorem by the total vector-meson 
nucleon cross section <tvn that will be specified below. The 
real part of the vector meson self energy SftlZy (M) in the medium 
might be calculated by dispersion relations, however, QCD sum 
rules for cc vector states suggest that the 'mass shift' 5R7T/(2My) 
at normal nuclear matter density is only in the order of a few 
MeV [ |l9| due to a small coupling of the c, c quarks to the nu- 
clear medium. We thus will assume SRTTy = throughout the 
following study. 



2.1 Vacuum properties 

Whereas the pole masses of the vector mesons My are approx- 
imately known e.g. from e + e~ annihilation [ |20| ] this does not 
hold for the total widths and especially branching ratios. To this 
end we have to model the vacuum decays, that enter the spec- 
tral function ([!]) via ([?]), using available information from the 
PDG pl"|]. In each channel c the relative decay to mesons mi 
and 77i2 is described by a matrix element (squared) and phase 
space, i.e. 

r c ~ \M c \ 2 P l (5) 
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that provide a reasonable description to the e + e formation 
data from Ref. [J2"o|], The resulting spectral function S(M) = 
2MA(M) for the #(4040) (with a pole mass of 4.08 GeV) is 
shown in Fig. 1 (upper part) in terms of the thick solid line 
as well as the mass dependent branching ratios for the differ- 
ent channels c (including the phase-space factor ~ p^). It is 
clearly seen that for the ratios (||) the D*D* decay opens up at 
M w 4.03 GeV and dominates the spectral function for higher 
masses. Furthermore, the decays to D S ,D* are strongly sup- 
pressed such that this resonance does not qualify for the pro- 
duction of (sc) or (cs) mesons. 

Apart from the total width of 80 ±25 MeV the properties of 
the #(4160) are not well known experimentally. For the decay 
channels we thus adopt a 'counting rule' for high energy open 
charm production that results from PYTHIA calculations (cf. 
Fig. 1 in Ref. [g), i.e. 
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(7) 



This is presently nothing but an educated guess (at low ener- 
gies) and must be controlled by future data. The resulting spec- 
tral function S(M) = 2MA(M) for the #(4160) is shown in 
Fig. 1 (lower part) in terms of the thick solid line as well as the 
mass dependent branching ratios for the different channels c. It 
is clearly seen that for the ratios (0) the D*D* decay still dom- 
inates the spectral function, however, the decays to D s , D* are 
no longer so strongly suppressed as in case of the #(4040). 

The production of the vector mesons in pp then can be de- 
scribed by Breit-Wigner resonance formation on the basis of 
the spectral function A(M 2 ) ([!]) - employing the proper kine- 
matics - provided that the branching of pp — > V is known, 



ay(M) = -Br pP ^v4MA v (M)r v ot {M) 
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Fig. 1. The spectral functions S(M) for the #(4040) and ^(4160) 
(thick lines, not normalized) and the mass differential branching ratios 
for both vector mesons within the model described in the text. 



where the factor 3/4 stems from the ratio of spin factors and 
k denotes the momentum of the p (or p) in the cms. In fact, 
the branching ratio is rather uncertain for the cases of interest 
here. Following our previous work [ |r7[ | we adopt Br(pp — > 
#(3770)) w 2 x 10~ 4 which is very similar to the branching 
of the #(2S) state. Due to a lack of detailed knowledge we 
assume this branching ratio also to hold for the #(4040) and 
#(4160). The resulting cross sections for the mesons V in pp 
reactions are displayed in Fig. 2 (upper part) as a function of 
the antiproton kinetic energy Tp in the laboratory reflecting the 
resonance structure dlb and the kinematics from m). 



2.2 pA reactions 

In order to simulate events for the reaction pA — > V we use the 
Multiple Scattering Monte Carlo (MSMC) approach. An ear- 
lier version of this approach - denoted as Intra-Nuclear Cas- 
cade (INC) model - has been applied to the analysis of -q and oj 
production in pA and pA interactions in Refs. [23 24| , p5| ,f26||. 
The production of the hidden charm vector resonances on nu- 
clei then can be evaluated within the MSMC by propagating the 
antiproton in the nucleus and folding the elementary pp — ► V 
production cross section with the local momentum distribution 




9 10 

T(p) [GeV] 

Fig. 2. The calculated cross sections for #(3770), #(4040) and 
$•(4160) for pp (upper part) and p + 12 C reactions (lower part) as 
a function of the antiproton kinetic energy in the laboratory. 



at the annihilation point. We have calculated the latter quantity 
in the local Thomas-Fermi approximation. 

The numerical results for the formation cross sections of 
the resonances V on 12 C are shown in the lower part of Fig. 
2 as a function of Tp indicating that the resonance formation 
is further smeared out by Fermi motion and the maximum in 
the cross section becomes reduced inspite of the larger p an- 
nihilation cross section on 12 C. The suppression factor in the 
maximum of the cross section - per elementary reaction - is of 
the order (see Ref. Q) 



S.F. 



Trr R m N 
(k F M R ) • 



(9) 



where Fr, Mr denote the vacuum width and mass of the pro- 
duced meson, whereas hp is the target Fermi momentum. This 
leads to factors of - 0.075, 0.14 and 0.21 for #(3770), # (4040) 
and *^ (4160) per pp reaction, respectively. On the other hand, 
the production on nuclei roughly scales with the target proton 
number as Z o e (see below) such that the cross section for 12 C 
in the maximum for the #(4160) is roughly the same as in 
pp annihilation, while it is lower for the narrower resonances 
#•(3770) and #(4040). 

Necessary parameters for a Monte Carlo (MC) simulation 
of rescattering are the elastic and inelastic VN scattering cross 
sections and slope parameters b for the differential cross sec- 
tions da e i jdt, which are approximated by 



d(J e i 
~~dt 



exp(M), 



(10) 



where t is the momentum transfer squared. These parameters 
as well as the masses of the rescattered particles determine the 
momentum and angular distributions of the particles in the fi- 
nal state. As in our previous study [ |l7[ | we use b= 1 GeV~ 2 
for D,D + N and VN scattering as an educated guess. Fur- 
thermore, the inelastic cross sections of the vector mesons with 
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nucleons have to be specified. Since the relative momenta in the 
V — N system for rescattering are in the order of a few GeV 
one might use high energy geometric cross sections for the dis- 
sociation to open charm mesons. For all the three resonances 
studied here we adopt an inelastic cross section of 1 1 mb as 
well as a e i = 1 mb. The inelastic cross section of 1 1 mb is 
larger by about a factor of 2 in comparison to the J/#, #(25) 
cross sections (see introduction) due to the larger size of the 
higher lying resonances. 

As mentioned above, the resonances V are propagated with- 
out any mean-field potential. This approximation should be 
well fulfilled since the expected potentials are only in the or- 
der of a few MeV [ |l9| while the energy of vector mesons in 
the lab. frame is a couple of GeV. Furthermore, due to the high 
mass of the resonances we neglect any intrinsic formation time 
tf for the resonances. 

The resonances are propagated with their actual mass M 
- selected by MC according to the spectral function (Q) - with 
velocity pv / Ey and decay in time according to the differential 
equation 

1 ' (ID 



mentally from K + p and K~p reactions. Whereas the K + (us) 
scatters only elastically with nucleons at low momentum - since 
the s cannot be exchanged with a light quark of a nucleon - the 
K~(su) cross section is dominated by resonant s-quark trans- 
fer leading to strange baryons such as A or £ hyperons. Similar 
relations are expected to hold for the D, D mesons, where es- 
pecially the light quark contribution should give much larger 
cross section on nucleons than the cc vector resonances. This 
analogy is based on SU(4)fi avor symmetry [ |3l"| ] which, how- 
ever, might be broken substantially in view of the different ge- 
ometrical sizes of K- and D-mesons. At present this is an open 
question, which has to be settled by experiment. 

For our calculations we adopt the following cross sections 
- taken as constants in the momentum regime of interest - 



dP v 
dt 



-~rZ t (M)P v (t) 

7 



with the total width (||), while 7 denotes the Lorentz 7-factor. 
Their decay to DD, DD* ,D*D,D*D*, D S D S , D*D„ ,D S D* S , 
D*D* is, furthermore, selected by Monte-Carlo according to 
the mass differential branching ratios displayed in Fig. 1. 

Due to the rather moderate collision rates involved (for Oy l N 
w 12m&) the life time of the vector mesons #(3770), #(4040) 
and #(4160) is not reduced sizeably in the nuclear medium 
relative to the vacuum (assuming the open charm mesons not 
to change their properties in the medium). This reduction then 
amounts to about 15% for the #(3770) , to - 8% for the #(4040) 
and to ~ 5 % for the #(4160), respectively. On the other hand, 
the D, D mesons might change their effective mass in the nu- 



clear medium as advocated in Refs. [28 , 29LBQ]. Especially for 



dropping D, D effective masses in the medium the decay width 
of the vector mesons becomes enhanced at finite density pa- 
This effect is most pronounced for the #(3770) since its pole 
mass is close to the D, D threshold of ~ 3.73 GeV. A drop of 
the D, D meson masses both by 50 MeV at density po would 
imply an increase of the decay width by about a factor of 2.5 ac- 
cording to phase space. Thus in-medium modifications of the 
open charm mesons might alter the life time of the #(3770) 
more effectively than collisional broadening. However, in- 
medium potentials are most pronounced at low momenta and 
decrease in size with increasing momentum. Thus it is very 
questionable if D, D mesons will have sizeable potentials for 
relative momenta of 3-5 GeV/c with respect to the nucleus at 
rest. Thus, in the following, we do not speculate any further 
on such medium mass modifications of D, D mesons at high 
relative momentum and continue our calculations for vacuum 
properties of the open charm system. 

The open charm mesons from the decaying vector mesons 
have a light quark q = (u, d) or strange s quark content apart 
from the c or c quark. In the constituent quark model we get 
D+ = (cd), D~ = (cd), D° = (cu), D° = (cu), D+ = (cs), 
D~ = (cs) and the same composition for the related vector 
states. Now light quark exchanges with nucleons ((uud) or 
(udd)) have different strength as it is well established experi- 
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(12) 



Here the inelastic cross sections of D-mesons refer to c-quark 
exchange reactions with nucleons. Furthermore, charge exchange 
reactions like D + n <-> D°p or D~p «-> D°n are incorpo- 
rated with a constant cross section a qexc . = 12 mb. For our 
exploratory study we assume the same cross sections for the 
related open charm vector mesons D* ,D* etc. 



3 Resonance production in pA reactions 

Within the model described in Section 2 we are now in the po- 
sition to calculate the cross section of the cc resonances in pA 
reactions for all targets of interest and can obtain the informa- 
tion about in-medium resonance decays (for densities pa > 
0.03/m~ 3 ) or vacuum decays, respectively. In this respect we 
show in Fig. 3 the calculated cross section for the #(3770) at 
T iab = 5.7 GeV for the targets 12 C, 64 Cu, 108 Ag and 207 Pb as 
a function of target mass A. The solid line represents the power 
law 

o* - A a (13) 

with a = 0.55, which nicely describes the total #(3770) yield. 
The coefficient a is less than 2/3 - as expected from the total p 
annihilation cross section - due to fact that only reactions with 
protons can lead to uncharged resonances and the Fermi smear- 
ing is more effective in heavy than in light nuclei. The fraction 
of #(3770), that decay to DD inside the nucleus (D ins) is 
only a few % in case of the 12 C target, however, reaches about 
50% for the Pb target. Thus for a heavy target the in-medium 
properties of the #(3770) may well be studied experimentally 
by dilepton spectroscopy exploiting the dilepton branching of 
7 x 10~ 4 of this resonance. However, as mentioned above, a 
collisional broadening of its spectral function by ~ 15 % is not 
expected to provide a clear signal relative to its vacuum decay 
unless the experimental mass resolution is in the order of 1-2 
MeV. 

Apart from the resonant production of open charm mesons 
they may also be produced directly as (cq) and (cq) pairs in 
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Fig. 3. The calculated cross sections for <?(3770) at T lab = 5.7 GeV as 
a function of the target mass A. The dotted line (D outs) denotes the 
contribution from 'outside' resonance decays whereas the dash-dotted 
line (D ins) shows the resonance contribution from 'inside' decays. 
The solid line reflects the scaling law ( pj| ) for a = 0.55 for the to- 
tal contribution from the resonance decay. The dashed line (D nonr) 
corresponds to the estimated nonresonant D, D cross section. 



Fig. 4. The calculated cross sections for <?(4040) at T iab = 6.8 GeV as 
a function of the target mass A. The dotted line (D outs) denotes the 
contribution from 'outside' resonance decays whereas the dash-dotted 
line (D ins) shows the resonance contribution from 'inside' decays. 
The solid line reflects the sum of 'outside' and 'inside' resonance de- 
cay contributions. The dashed line (D nonr) corresponds to the esti- 
mated nonresonant D, D cross section. The lower two lines give the 
'outside' and 'inside' decay to D s ,D a , respectively. 



pN annihilation. This channel strongly dominates in annihila- 
tion reactions on neutrons in the nucleus since the charmonium 
production is forbidden by charge conservation at low energy 
above threshold. At higher invariant energies a pion might bal- 
ance the charge in the ty+w final channel. For our estimates we 
have employed the Regge-model analysis of Ref. [ f32| ] for the 
elementary pN — > D,D cross section as also used by Sibirt- 
sev et al. in Refs. [ ^8[p9[ |. The resulting cross section of non- 
resonant open charm production is shown in Fig. 3 as a func- 
tion of the target mass A in terms of the dashed line (D nonr) 
and demonstrates that the nonresonant cross section is about 
the same at this energy as the resonant production channel via 
5^(3770) production and decay. Thus, in-medium modifications 
of the ^(3770) will be hard to detect experimentally via the 
DD invariant mass spectrum (see below). 

In analogy to Fig. 3 we show in Figs. 4 and 5 the calculated 
cross section for the <f(4040) at T iab = 6.8 GeV and ^(4160) 
at Ti ab = 7.7 GeV for the targets 12 C, 64 Cu, 108 Ag and 207 Pb 
as a function of target mass A. Due to the shorter life times of 
these resonances the fraction of in-medium resonance decays 
increases up to 65 % and 75%, respectively, for a P6-target. 
However, for the rescattering cross sections involved (cf. Sec- 
tion 2) the life times in the medium are shortened only by 8% 
and 5%, respectively. This very moderate change of the spectral 
functions due to in-medium interactions will be hard to see by 
dilepton spectroscopy except for setups with excellent invariant 
mass resolution. 



For the ^(4040) the nonresonant D, D cross section (upper 
dashed line; D nonr) becomes larger than the resonant cross 
section (solid line) due to the higher invariant energy yfs in- 
volved in the elementary reactions at Ti ao = 6.8 GeV. Fur- 
thermore, the resonant decays to D s , D s are suppressed by al- 
most 3 orders of magnitude in case of the !f - (4040) which is 
not favorable for studying D s , D s propagation in the nuclear 
medium. For the 1^(4160) an even larger nonresonant contribu- 
tion is expected which we do not attempt to estimate explicitly 
since no educated guess for the additional channels DDtt + X 
is available. However, the perspectives for D s , D s dynamics 
become more promising due to larger branching ratios which 
are a consequence of the higher invariant energy above the 
threshold and a more favorable matrix element. Especially for 
the heaviest targets almost 50% of the D s , D s mesons from 
<F (4160) decay appear inside the nucleus such that rescattering 
effects can be explored. Such rescattering phenomena can also 
be tested for nonresonant D S ,D S production channels; how- 
ever, their cross section is (at present) hard to estimate as men- 
tioned before. 



4 Open charm propagation and rescattering 

The open charm mesons produced by resonance decay or non- 
resonant channels may rescatter in the nuclear medium either 
elastically or inelastically. In case of inelastic reaction channels 
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Fig. 5. The calculated cross sections for *'(4160) at T lab = 7.7 GeV as 
a function of the target mass A. The dotted line (D outs) denotes the 
contribution from 'outside' resonance decays whereas the dash-dotted 
line (D ins) shows the resonance contribution from 'inside' decays. 
The solid line reflects the sum of 'outside' and 'inside' resonance de- 
cay contributions. The lower two lines give the 'outside' and 'inside' 
decay to D 3 ,D 3 , respectively. 




P(D),P(D) [GeV/c] 

Fig. 6. The calculated momentum distributions of D, D mesons for 
p + Pb reactions at Ti a b = 5.7 GeV. The hatched area reflects open 
charm mesons from 1^(3770) decays in vacuum, the dotted (Dr) and 
dash-dotted (Dr) lines correspond to resonance decays in the target nu- 
cleus, respectively. The solid and dashed lines display the total spectra 
for D and D mesons including the nonresonant production channels 
(see text). In the calculations elastic and inelastic scatterings of the 
$■(3770) and D, D mesons have been taken into account. 



there may be charge exchange reactions such as D + n <-> D°p 
or D~p <-> D°n or charm flavor exchange reactions DN — > 
A c tt or S c tt in analogy to the strangeness exchange reactions 
K N — > Air or Stt. Whereas charge exchange reactions - in- 
duced by (u, d) quark exchange - do not change the ratio D/D 
in nuclei, the charm exchange reaction leads to a loss of c 
quarks in the mesonic sector and thus to a decrease of the D/D 
ratio with the target mass number A. 

In order to obtain some idea about the momenta of open 
charm mesons in the laboratory we display in Fig. 6 the mo- 
mentum distribution of D and D mesons - after full rescattering 
- for p+Pb at Ti a b = 5.7 GeV. At this energy the resonance con- 
tribution stems from the >P(3770); its 'outside' decay contribu- 
tion is shown in Fig. 6 by the hatched area, which extends from 
2.2 GeV/c to about 4.2 GeV/c. The 'inside' resonance contri- 
bution (denoted by Dr and Dr, respectively) is seen to be broad- 
ened significantly by elastic scattering and charge exchange re- 
actions showing also a sizeable net absorption of D mesons in 
case of the Pb target of ~ 30%. On the other hand, the D, D 
mesons from nonresonant production channels signal larger in- 
medium rescattering effects which is seen by comparison to 
the total D (solid line) and D (dashed line) momentum spec- 
tra. The latter phenomenon is due to the fact that rescattering 
of open charm mesons in case of nonresonant production chan- 
nels may proceed earlier than in case of the ^(3770) — > DD 
channel since the 5^(3770) on average is propagating for a cou- 
ple of fm/c before its decay. 



In order to quantify the observable consequences from D, D 
rescattering and absorption we have made a cut on the labora- 
tory momenta Po,Pg < 2.2 GeV/c to exclude the vacuum 
decays of the ^(3770) and to reduce the fraction of unscat- 
tered D, D events. The calculated cross sections for resonant 
(l.h.s.) and nonresonant (r.h.s.) production channels are shown 
in Fig. 7 (upper part) for D and D mesons, respectively, as 
a function of the target mass A at Ti ab = 5.7 GeV. Even for 
the low momentum cut the cross sections are in the order of 
50 - 100 nb for heavy targets, which should be feasible at the 
HESR with a p luminosity of — 2x 10 32 . The resulting ratio 
D to D mesons is displayed in the lower part of Fig. 7 for res- 
onant (l.h.s.) and nonresonant (r.h.s.) productions separately. 
The constant line (=1) results trivially when including charge 
exchange reactions, however, discarding c-quark exchange re- 
actions. As noted above, the D-meson absorption is larger in 
the nonresonant case and may reach 40% for the P6-target 
within the cross sections specified in Sec. 2. Vice versa, exper- 
imental data on this ratio (employing the same low momentum 
cut) as a function of mass A will allow to set stringent bounds 
on the size of the c-quark exchange cross section in the nuclear 
medium. 

Without explicit representation we note that the perspec- 
tives of measuring D, D rescattering in heavy nuclei from the 
resonance decay of the !f - (4040) are similar to the previous case 
discussed above and also close to the results for the 5^(4160) 
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Fig. 7. Upper part: The calculated cross sections for D and D mesons 
for p + A reactions at Ti a b=5.7 GeV including a low momentum 
cut Pd, Pjj < 2.2 GeV/c. The l.h.s. shows the contribution from 
in-medium resonance decays whereas the r.h.s. corresponds to non- 
resonant production channels. Lower part: The ratio D/D versus the 
target mass A for resonant (l.h.s.) and nonresonant (r.h.s.) production 
channels. The constant line (=1) results when neglecting c-quark ex- 
change channels. 
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(see below). Since the chances to measure additionally the re- 
scattering effects of D+ , D~ mesons are very low in view of 
the small partial decay width for the ^(4040) (cf. Fig. 1) we 
directly continue with a discussion of the ^(4160) formation 
and decay at Ti ao = 7.7 GeV. 

In this case (for p + Pb) the momentum distribution of the 
D, D, D+,D~ mesons becomes very broad and extends from 
~ to - 7 GeV/c (cf. Fig. 8) where the 'outside' ^(4160) 
decays (solid lines; D, D outside, Df, D~ outside) give open 
charm mesons essentially above <~ 2 GeV/c. Thus rescatter- 
ing events of the D-mesons are dominantly found for momenta 
less than 2 GeV/c which, unfortunately, are only a low frac- 
tion of the total open charm yield. Furthermore, the number 
of Df,D~ events in this low momentum range is down by 
more than 2 orders of magnitude relative to the D, D mesons. 
However, apart from the resonant production of charm meson 
pairs there will be an additional contribution from nonresonant 
D, D, , Dj mesons which will enhance the cross section 
also for low momenta, respectively. It is hard to provide reli- 
able estimates for the nonresonant contribution, but in view of 
Fig. 4 for Ti a b = 6.8 GeV this contribution might be higher by 
a factor 2-4 than the resonance contribution. 

As shown in Fig. 9 the cross sections for D, D mesons at 
Tiab = 7.7 GeV are above 50 nb for heavy targets (l.h.s, upper 
part) - even when including a low momentum cut of 2 GeV/c 
- showing a sizeable difference due to the absorption of D, D* 
mesons in the nucleus. The perspectives for D+,D~ mesons 



are less promising (r.h.s, upper part). Here only cross sections 
of <~ 0. 1 nb should be expected for a Pb target which might 
be enhanced due to nonresonant channels by a factor 2^1 (as 
estimated above). However, the relative absorption of D, D* 
mesons - as measured by the ratio to D, D* mesons - increases 
with target mass by up to 35% for A=207 (l.h.s., lower part) 
for the cross sections specified in Sect. 2. This decrease in the 
D/D ratio is roughly the same for the !f - (4160) decays as for 
the ^(3770) decays (Fig. 7) inspite of the shorter lifetime of 
the ^(4160). This is partly due to the different cut in momen- 
tum and the fact, that the average momentum of the produced 
^(4160) in the laboratory is higher. The ratio of D~ /Df ver- 
sus target mass A (r.h.s.; lower part) also shows a decrease due 
to the different cross sections for c and s quark exchanges with 
the nucleons of the target. The vertical lines in Fig. 9 show 
the statistical error bars in this ratio for the low momentum cut 
of 2 GeV/c that reflect the low number of D~ , Df mesons in 
this momentum range as obtained from the Monte Carlo sim- 
ulations. In principle, provided that sufficient statistics are ob- 
tained experimentally, it is possible to determine the relative 
strength of (u, d), s and c quark exchanges with nucleons from 
the mass dependence of the further ratios D/D+ and D/D~. 
We only point out this possibility without presenting explicit 
figures for these ratios since the statistics for Df is very small 
and the cross sections employed are rather uncertain so far. 
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Fig. 9. Upper part: The calculated cross sections for D and D mesons 
for p + A reactions at Ti a b=l '■! GeV including a low momentum cut 
Pd and Pg < 2.0 GeV/c. The l.h.s. shows the contribution from in- 
medium resonance decays for D, D whereas the r.h.s. corresponds 
to Df , D~ resonant production (multiplied by a factor 100). Lower 
part: The ratio D/D versus the target mass A (l.h.s.) and the ratio 
D~ /Df for resonant production. The constant line (=1) results when 
neglecting s and c quark exchange channels. 



5 Summary 

In this study we have explored the perspectives of measuring 
the dynamics of hidden charm vector mesons and open charm 
mesons in nuclei in antiproton induced reactions on nuclei. 
Such experimental studies will be possible at the high-energy 
antiproton storage ring HESR proposed for the future GSI fa- 
cility 0]. The Multiple Scattering Monte Carlo (MSMC) cal- 
culations performed have been based on the resonance produc- 
tion concept with resonance properties from the PDG or 
related estimates from high energy open charm branching ra- 
tios. We note, that these resonance properties are in line with 
#(3770), #(4040) and #(4160) production from e+e~ anni- 
hilation pc|, however, the detailed branching ratios for the 
#(4040) and #(4160) (cf. Fig. 1) are still quite uncertain. Thus 
our estimated uncertainty in the cross sections is in the order 
of a factor 2-3, which should be kept in mind when planning 
experiments along this line. We stress that the total branching 
ratios to pp should be measured first with high accuracy. De- 
viations from our calculations for hidden charm vector mesons 
in pp reactions than directly carry over to the cross sections 
on nuclei (Fig. 2). Nevertheless, even for a light nucleus such 
as 12 C cross sections from 70 - 150 nb should be expected at 
the proper resonance energies (Fig. 2). Our calculations sug- 
gest that the total yield of these vector mesons scales as ~ 
A0.55 w j tn taj-ggt mass A (Fig. 3). For a heavy nucleus like 
Pb the fraction of 'inside' decays to dileptons or D, D reaches 



- 50% for the #(3770), - 65 % for the #(4040) and - 75% 
for the #(3770), which is sufficient to address the question of 
in-medium properties of these resonances. However, for the 
rescattering cross sections adopted in Sect. 2, the in-medium 
modifications of their spectral functions is rather moderate, i.e. 
15%, 8% and 5% change in width, respectively, such that a 
mass resolution for dilepton pairs in the order of a few MeV 
will be necessary. This will be also the case for in-medium 
mass shifts of these vector mesons which are estimated to be 
in the order of a few MeV, only [ |l9[ ] . 

We have, furthermore, explored the perspectives of mea- 
suring the interactions of open charm mesons in nuclei that 
are produced by resonant vector meson decays or nonresonant 
channels, respectively. For the cross sections specified in (12) 
the effects of D, D rescattering and D-meson absorption in 
heavy nuclei can clearly be identified when imposing a low 
momentum cut in the laboratory of 2-2.2 GeV/c at Ti a \, = 5.7 
or 7.7 GeV, respectively. This low momentum cut excludes 
events with DD decays in vacuum or D, D mesons that did 
not rescatter in the target nucleus. The cross sections for such 
type of events with D or D mesons of low momenta are in 
the order of 50 nb for a heavy target like Pb. Such cross sec- 
tions are high enough to perform successful experiments with 
an antiproton luminosity of 2 x 1 32 as envisaged for the HESR 

. Experimentally, the ratio of D/D (at low momenta) as a 
function of target mass A allows to set stringent constraints 
on the Z?-meson absorption - or c-quark transfer - cross sec- 
tion at nuclear density po- Furthermore, at T; a fc w 7.7 GeV 
the production cross section of 7J+ , D~ mesons should be- 
come appreciable and be in the order of 20 nb for a Pb target 
when including only the #(4160) decays. Approximately 50% 
of these (es) and (cs) mesons should rescatter in a target of 
mass A sa 200 such that c and s-quark exchange reactions with 
nucleons can be measured by studying the mass dependence of 
Dj /D+ as well as Dj D~ and D/D/. at low momenta. How- 
ever, to exclude events from #(4160) vacuum decays as well 
as events with unscattered open charm mesons a severe cut on 
low momenta (< 2 GeV) should be performed. In this momen- 
tum range the cross section for D~ , Df drops to ~ 0. 1 nb even 
for a Pb target. It will thus be very hard to measure the various 
particle ratios as a function of mass A with good statistics. 
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